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Abstract 

We consider QCD radiative corrections to top-quark pair production at hadron 
colliders. We use the qx subtraction formalism to perform a fully-differential 
computation for this process. Our calculation is accurate up to the next-to-leading 
order in QCD perturbation theory and it includes all the flavour off-diagonal partonic 
channels at the next-to-next-to-leading order. We present a comparison of our 
numerical results with those obtained with the publicly available numerical programs 
MCFM and Top++. 
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The top quark (t) has a special role [I] in elementary particle physics. Being the heaviest 
known fundamental constituent, with a mass of about 173.3 GeV [2], it couples strongly to the 
Higgs boson and it is crucial to the hierarchy problem. Within the Standard Model (SM) the 
main source of top-quark events in collisions at hadron colliders is top-quark pair production. 
Many New Physics (NP) models predict the existence of top partners with masses close to the 
electroweak symmetry breaking scale, which exhibit similar properties as the top quark and can 
decay into it. Studying the production of tt pairs at hadron colliders can not only shed light on the 
nature of the electroweak-symmetry breaking but it also provides information on the backgrounds 
of many NP models. 

The theoretical efforts for obtaining precision predictions for top-quark pair production at 
hadron colliders started almost three decades ago with the calculation of the next-to-leading order 
(NLO) QCD corrections to the total cross section [3y5j and kinematical distributions [6] for this 
production process. The NLO calculations of the total cross section of Refs. were carried 
out numerically. The expressions in analytic form of the total partonic cross sectiorU at NLO 
were obtained in Ref. [TO] . Recently the calculation of the next-to-next-to-leading order (NNLO) 
QCD corrections to the tt total cross section was completed HH- Besides the total cross section, 
differential cross sections and more general kinematical distributions are of great importance for 
precision studies. For instance, the tt (forward-backward and charge) asymmetry has received 
much attention in recent years (see, e.g., Ref. |12j). The tt asymmetry, which is non vanishing 
starting from the NLO level [13], has recently been computed up to the NNLO level [Llj. Other 
NNLO results on differential distributions are starting to appear mm- 

This Letter is devoted to the NNLO (and NLO) QCD calculation of tt production. In partic¬ 
ular, we present the results of the first NNLO application of the q? subtraction formalism [T8] to 
the process of tt production in hadron collisions. 

At the partonic level, the NNLO calculation of tt production requires the evaluation of tree- 
level contributions with two additional unresolved partons in the final state, of one-loop contri¬ 
butions with one unresolved parton and of purely virtual contributions. The required tree-level 
and one-loop scattering amplitudes are known and they are the same amplitudes that enter the 
NLO calculation of tt + jet [19], the associated production of a tt pair and one jet. The purely 
virtual contributions depend on the two-loop scattering amplitudes and on the square of one- 
loop scattering amplitudes. The two-loop amplitude for tt production is partly known in analytic 
form [20] and its complete computation has been carried out numerically [21] . The square of 
one-loop scattering amplitudes is known [22] . 

The implementation of the various scattering amplitudes in a complete NNLO calculation at 
the fully differential (exclusive) level is a highly non-trivial task because of the presence of infrared 
(IR) divergences at intermediate stages of the calculation. In particular, these divergences do not 
permit a straightforward implementation of numerical techniques. Various methods have been 
proposed and used to overcome these difficulties at the NNLO level. The formalisms of antenna 
subtraction [2511211] and colourful subtraction E3EB] are more related to NNLO extensions of 
established NLO formulations [29l [3T| of the subtraction method. The Stripper formalism [321134] is 
a combination of the subtraction method with numerical techniques based on sector decomposition 
[351136]. Variants of the subtraction methods are the qx subtraction formalism [18] and the recently 

tA parametrization [7] of this analytic NLO result is implemented in the numerical programs Hathor [8] and 
TOP++ 0. 
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proposed N-jettiness subtraction |57H5^j . 


The NNLO computations in Refs. mm for tt production have been performed by using 
the Stripper method [32]. Parallcly, an ongoing effort is being carried out by using the antenna 
subtraction method [241140] , which led to the NNLO fully differential computation of tt production 
in the qq channel [TGlITT] at leading colour and including the light-quark contributions. 


The qx subtraction formalism [L8] is a method to handle and cancel the IR divergences at 
the NLO and NNLO level. The method has been successfully applied to the fully differential 
computation of NNLO QCD corrections to several hard-scattering processes [ I8H42H44] . The 
method uses IR subtraction counterterms that are constructed by considering and explicitly com¬ 
puting the transverse-momentum (qx) distribution of the produced final-state system in the limit 
qx —> 0. If the produced final-state system is composed of non-QCD (colourless) partons (e.g., 
leptons, vector bosons or Higgs bosons), the behaviour of the qx distribution in the limit qx —i > 0 
has a universal (process-independent) structure that is explicitly known up to the NNLO level 
through the formalism of transverse-momentum resummation [45] . These results on transverse- 
momentum resummation are sufficient to fully specify the qx subtraction formalism for this entire 
class of processes. Therefore, up to now, the applications of the qx subtraction formalism have 
been limited to the production of colourless high-mass systems in hadron collisions. In this Letter 
we present first results on the application of the qx subtraction method to the NNLO computation 
of heavy-quark production in hadron collisions. To this purpose, we use the recent progress on 
transverse-momentum resummation for heavy-quark production [461448] . We exploit the formula¬ 
tion of transverse-momentum resummation in Ref. [48] that includes the complete dependence on 
the kinematics of the heavy-quark pair. This dependence and, in particular, the complete con¬ 
trol on the heavy-quark azimuthal correlations are essential (see below) to extract all the NNLO 
counterterms of the qx subtraction method. Although the structure of transverse-momentum 
resummation for heavy-quark production is fully worked out up to the NNLO level, the explicit 
NNLO results for the hard-virtual factors [48] in the flavour diagonal partonic channels qq —> tt+X 
and gg ^ tt + X (X denotes the unobserved inclusive final state) are not yet known. Therefore, 
in the NNLO calculation of this paper we present numerical results for all the flavour off-diagonal 
channels ab tt + X, with ab = qg(qg),qq(qq) J qq , (qq'),qq , (qq') ( q and q' denote quarks with 
different flavour). 


The differential cross section da tt for the inclusive production process pp(pp) —$■ tt + X is 
computable by convoluting the corresponding partonic cross sections dd^ b of the various partonic 
channels with the parton distribution functions (PDFs) of the colliding hadrons pp (pp). According 
to the q T subtraction method [IB], the (N)NLO partonic cross section can be written as 


da 


tt 

(N)NLO 


n/tt 

n (N)NLO 


® d&LO 


+ 


d ^N)L 0 ~ dd 


tt, CT 
( N)NLO 1 


( 1 ) 


where da^yto the R+jet cross section at (N)LO accuracy. Applying Eq. ([!]) at NLO, the 
leading-order (LO) cross section da ^Q Jct can be directly obtained by integrating the corresponding 
tree-level scattering amplitudes. Applying Eq. (pQ) at NNLO, da^Q can be evaluated by using any 
available NLO method (e.g., Refs. [29H3T] 1 to handle and cancel the corresponding IR divergences. 
Therefore, dcr^pLo is IR finite provided qx ^ 0. 


The square bracket term of Eq. (JT]) is IR finite in the limit qx —> 0, but its individual contribu¬ 
tions, da^^Q and da are separately divergent. The IR subtraction counterterm da^'^ LO 
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is obtained from the (N)NLO perturbative expansion (see, e.g., Ref. [491I5U] ) of the resummation 
formula of the logarithmically-enhanced contributions to the qx distribution of the tt pair [06H48]: 
the explicit form of ckpZ)p LO can ^* e coiaipletely worked out up to NNLO accuracy. For example, 
at the NLO, the explicit expression of da N ’ LO in the partonic channel ab —$■ tt + X is 


, - tt, ct _ Q-s 

aa NLO ab~ n 

c=q,q,g 


yR) I y-'(l)**— new> \ 
cc\ — ab cd — cib J 


dqj, 

JP 


<g> da tt 


LO cc 


( 2 ) 


where as is the QCD coupling, M is the invariant mass of the produced tt pair and da^ 0ab is the 
LO partonic cross section. The expression ([2]) involves convolutions (which are denoted by the 
symbol with respect to the longitudinal-momentum fractions Z\ and Z 2 of the colliding partons 
c and c in dcr/ 0 c5 . The integration variable q T in Eq. ([2]) corresponds, in the limit q T —* 0, to the 
transverse momentum of the produced tt pair in the cross section da^p }et on the right-hand side 
of Eq. (pQ). The function E^_ ob enters into the qx subtraction method [18] for hard-scattering 
production of a generic final-state system. Its explicit form is [191150] 


s S-o6(*i»*2;gr/A0 


- ^A^dcaScbdi 1 - Zi)5(l - z 2 ) h{qr/M) - 
+ S ca 6( 1 - Zl ) P£\z 2 ) + 63,6(1 - z 2 ) P£>(z 1 ) 


S ca 3 cbd(l - zi) 5 (l - z 2 ) B ^ 
h(qr/M) , (3) 


and it derives from the small-g^ singular behavior of the qx cross section for the production of a 
colourless system in the partonic cc production channel. The coefficients AP and B P are the first- 
order resummation coefficients for transverse-momentum resummation (AP = Cf, Ap = Ca, 
Bp = —3/2 Cf , Bp = —(11/6 Ca -%/3)). The functions Pp(z) are the lowest-order DGLAP 
kernels (the overall normalizazion is specified according to the notation in Eq. (41) of Ref. [49]). 
The functions Ipqx/M ) (k = 1,2), which appear in Eq. Q, encapsulate the singular behavior at 
small qx, and they are explicitly given in Appendix B of Ref. [12]. The other function E^^ afe new 
in the round-bracket factor of Eq. ([2]) is due to soft radiation and it is an additional term that is 
specific of the qx subtraction method for the case of heavy-quark pair production. This function 
reads 


, - (Mct-tf I (r, (1) + rg )f ) 

1 , 22 ; q T /M) = -Uf( 1 - Zi)«(i - 22 )- —P - Uqr/M ), 

(4) 

where r) 1 ^ is the first-order term of the soft anomalous dimension for transverse-momentum re¬ 
summation in heavy-quark production and its explicit expression is given in Eq. (33) of Ref. [48j |. 
This soft anomalous dimension is a colour space matrix that acts onto the colour indices of the four 
partons {c, c, t, t} in the Born level scattering amplitude \M. C c^it) of the partonic process cc —* tt. 
The colour space notation is specified in Ref. [48] and, in particular, \A4 c5 ^. t p = (Xi c5 ->tt\M.cc-+tt) 
denotes the colour summed square amplitude that contributes to dct/ 0cg , whereas the factor 
(M^ t i\(rP + TP^lMcc^ti) embodies colour correlation terms with a definite kinematical de¬ 
pendence. 

The first-order hard-collinear (IR finite) counterterms are also completely known [41)1148] 

for all the partonic channels. The second-order (IR finite) counterterms 'Hnnlo are n °f y e f 
fully known. However, 'H% NLO can be explicitly determined for all the flavour off-diagonal par- 
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tonic channels. In these off-diagonal channels, H% NLO embodies process-dependent and process- 
independent contributions. The process-dependent contributions to derive from the knowl¬ 

edge of the one-loop virtual amplitudes of the partonic processes qq —¥ ti and gg —> ti, and from 
the explicit results on the NLO azimuthal correlation terms in the transverse-momentum resum¬ 
mation formalism [48] (see, in particular, Eq. (25) in Ref. |@8j and accompanying comments). 
The process-independent contributions to are analogous to those that contribute to Higgs 

boson [18] and vector boson [43] production, and they are explicitly known [[45ll51Ii52 j. 

Having discussed the content of Eq. ©, we are in a position to apply it to tt production 
and to obtain the complete NLO results plus the NNLO corrections in all the flavour off-diagonal 
partonic channels. Our NLO implementation of the calculation has the main purpose of illustrating 
the applicability of the qx subtraction method to heavy-quark production and, in particular, of 
cross-checking the qx subtraction methodology by numerical comparisons with NLO calculations 
performed by using more established NLO methods. Our NNLO results on ti production represent 
a first step (due to the missing flavour diagonal partonic channels) towards the complete NNLO 
calculation for this production process. Up to NLO accuracy our numerical implementation is 
based on the scattering amplitudes and phase space generation of the MCFM program [53j , 
suitably modified for qx subtraction along the lines of the corresponding numerical programs for 
Higgs boson [18] and vector boson [43] production. At NNLO accuracy the tt+jet cross section 
is evaluated by using the Munich code [53], which provides a fully automated implementation 
of the NLO dipole subtraction formalism [30ll3l] as well as an interface to the one-loop generator 
OpenLoops [55] to obtain all the required (spin- and colour-correlated) tree-level and one-loop 
amplitudes. For the evaluation of tensor integrals we rely on the COLLIER library [56], which 
is based on the Denner-Dittmaier reduction techniques of tensor integrals and on the scalar 
integrals of Ref. [58]. In OpenLoops problematic phase space points are addressed with a rescue 
system that uses the quadruple-precision implementation of the OPP method in CutTools [59] 
with scalar integrals from OneLOop [60] . 




Figure 1: The invariant mass (left) and rapidity (right) distributions of the ti pair at the LHC 
(i/s = 8 TeV) computed at NLO accuracy. Comparison of our results (blu) with the MCFM 
results (red). The lower panel presents the ratio of our results over the MCFM results. 
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Figure 2: The rapidity (left) and transverse-momentum (right) distributions of the top quark at 
the LffC (y/s = 8 TeV) computed at NLO accuracy. Comparison of our results (blu) with the 
MCFM results (red). The lower panel presents the ratio of our results over the MCFM results. 


We start the presentation of our results by considering pp collisions at y/s = 8 TeV. We use 
the MSTW2008 [6Tj PDFs with the QCD running coupling as evaluated at each corresponding 
order (i.e., we use (n + l)-loop as a t N”LO, with n = 1,2). The pole mass of the top quark is 
m t = 173.3 GeV. The renormalization and factorization scales, pp and pf, are fixed at prt = Pf = 
m t . 


In Figs. H] and [2] we compare the NLO differential distributions obtained by using MCFM 
(which implements the dipole subtraction method [3Q)|3T]) with those obtained by using our 
numerical program^. In particular in Fig. |Tj we consider the invariant mass ( m t t ) distribution 
(left) and the rapidity ( y t j) distribution (right) of the tt pair. In Fig. [2] we consider the rapidity 
{yt) distribution (left) and the transverse-momentum {pr,t) distribution (right) of the top quark. 
We clearly see that the distributions obtained with qp subtraction are in excellent agreement 
with those obtained with MCFM. We have checked that the agreement persists also for different 
choices of pp and pp- 

We now move to consider the NNLO contributions and, in particular, we compute the total 
cross section for tt production. In Table |Tj we report our results and we compare them with 
the corresponding results from the numerical program TOP++ [9j, which implements the NNLO 
calculation of Ref. HU • Specifically, we report the complete NLO results and the 0(a|) contri¬ 
butions to the NNLO cross section from the flavour off-diagonal partonic channels ab —» tt + X. 
The contribution from all the channels with ab = qg , qg is labelled by the subscript qg , and the 
contribution from all the channels with ab = qq , qq , qq\ qq', qq ', qq' is labelled by the subscript 
q{q)q'. From Table [T] we see that the results obtained by using qp subtraction are in agreement 
with those from TOP+ + . We note that the numerical uncertainties of our 0(a g) results are at 

*These NLO results are obtained by using MCFM-V7.0 with about 9 hours of run on an Intel Xeon 2.4 Ghz, 
corresponding to a total of 2 * 10 6 integral evaluations. Our NLO implementation of qp subtraction requires about 
30 times higher statistics to obtain comparable results and the corresponding runtime is a factor of 3 larger. 
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the 2% level. This is due to the fact that in both the qg and q(q)q' channels at y/s = 8 TeV there 
is a strong quantitative cancellation between the contributions of the two terms in the right-hand 
side of Eq. ([T]) (the term that is proportional to an d fl ie t erm i n the square bracket). The 

numerical uncertainties of our 0{a§) calculation can be reduced by considering different centre- 
of-mass energies. In particular, the numerical cancellation that we have mentioned is reduced by 
decreasing the centre-of-mass energy of the collision. We have computed the total cross section 
for pp collisions at y/s = 2 TeV and we report the corresponding results in Table [2j We note 
that the numerical agreement between our calculation and the TOPH—K result is still satisfac¬ 
tory, and the numerical uncertainties of our 0(a|) results are reduced below the 1% level. Using 
TopH —b we can compute the complete NNLO result and we note that the flavour off-diagonal 
partonic channels contribute to about 10% of the total result at 0(a g) for both collision energies 
considered in Tables Q] and [2] 


Cross section [pb] 

NLO 

0(a$)\ qg 

0(<4) 9 ( 9 ) 9 ' 

qx subtraction 

226.2(1) 

-2.25(5) 

0.151(3) 

TOP++ 

226.3 

-2.253 

0.148 


Table 1: Total cross sections for tt production. NLO and (partial) NNLO results from qx subtrac¬ 
tion compared with the corresponding results from TOPH—|- for pp collisions at y/s = 8 TeV. 


Cross section [fb] 

NLO 

o(<4) I 9.9 

°( a s)\q(qW 

qx subtraction 

7083(3) 

-61.5(5) 

1.33(1) 

Top++ 

7086 

-61.53 

1.33 


Table 2: Total cross sections for tt production. NLO and (partial) NNLO results from qx subtrac¬ 
tion compared with the corresponding results from TOPH —\- for pp collisions at y/s = 2 TeV. 

We have presented the first application of the qx subtraction method to top-quark pair produc¬ 
tion at hadron colliders. Our implementation is based on the formulation of transverse-momentum 
resummation for heavy-quark production of Ref. [38], which includes the complete dependence on 
the kinematics of the heavy-quark pair. Our computation is accurate at NLO in QCD perturba¬ 
tion theory and it includes all the off-diagonal partonic channels at NNLO accuracy. At NLO we 
have compared our results for various kinematical distributions with those obtained by using the 
MCFM program, and we find good agreement. At NNLO our results for the tt total cross section 
agree with the corresponding results obtained by using the TOP++ program. The extension of 
our NNLO computation to include the missing qq —> tt + X and gg —> tt + X channels requires 
the evaluation of the second-order hard-collinear functions 'H^lo 03], and an implementation 
of the two-loop virtual amplitudes, which, at present, are known only in numerical form [21 j. 

The computation that we have performed in this paper can straightforwardly be extended to 
the production of massive-quark pairs of different flavour (e.g. bottom-quark pair). The extension 
of the method to production processes with massless coloured particles in the final state (e.g. 
inclusive dijet production) is definitely non trivial and it would require additional theoretical 
advancements. 
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